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Abstract 

Wc present a numerical study of the magnetic structure of nanostruc- 
tured iron fluoride, using the Monte- Carlo-Metropolis simulated annealing 
technique and a classical Heisenberg Hamiltonian with a superexchange an- 
gle dependence. The parameters are adjusted on experimental results, and 
the atomic structure and topology taken from a previous atomistic model 
of grain boundaries in the same system. We find perfect antiferromagnetic 
crystalline grains and a disordered magnetic configuration (speromagnetic 
like) at the grain boundary, in agreement with experimental findings. Both 
the lowest magnetic energy and the rate of magnetic frustration are found to 
be dependent on the relative disorientation of crystalline grains, i.e. on the 
cationic topology. By simulating hysteresis loops, we find that the magneti- 
zation rotation is not spatially uniform. We conclude on possible extensions 
of the model. 

Keywords: nanostructuration, iron fluoride, grain boundaries. Metropolis, 
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1. Introduction 

During the last decade, strong efforts have been devoted studying the 
magnetic properties of nanostructures. This growing interest is due to their 
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unusual and tunable physical properties, which are strongly influenced by the 
effect of conflnement resulting in a large contribution of surface and interface 
effects [1]. Understanding the correlation between magnetic properties and 
nanostructure involves collaborative efforts between chemists, physicists and 
materials scientists to study both fundamental properties and potential ap- 
plications [2]. This can be done experimentally using diffraction techniques, 
local probe techniques such as Mossbauer spectrometry, nuclear magnetic res- 
onance, electron microscopies and magnetic measurements. In addition, one 
can use numerical modeling techniques based on ab initio or phenomenolog- 
ical calculations such as the classical Monte Carlo methods with the Heisen- 
berg model, as in [31 IH El E] for instance. The key issue for understanding 
magnetic macroscopic properties, such as magnetization or susceptibility, 
would be to investigate the contributions arising from the nanograins and 
the role played by the nanograin surface and the interface between grains 
generally named grain boundary (GB). GBs in magnetic nanostructures fea- 
ture a disordered atomic structure and a spin-glass-like behavior and have 
a chemical composition which might differ strongly from that of the corre- 
sponding nanograin [3 El |9] . 

To investigate the magnetic behavior of nanostructures and the mutual in- 
fluence of the grains and grain boundaries, we selected ferric fluoride (FeF^,), 
because this ionic phase displays several advantages. Indeed, it exhibits 
a fascinating polymorphism with three different crystalline phases and two 
amorphous varieties, the structures of which result from ordered and disor- 
dered packing of corner-sharing octahedral units respectively. Those rather 
different cationic topologies originate coUinear and non-coUinear magnetic 
arrangements in conjunction with the antiferromagnetic nature of the su- 
perexchange interactions. Iron fluoride is therefore considered as an excel- 
lent simple model system to illustrate the concept of topological magnetic 
frustration. The more stable crystalline form of FeF^, which is the rhombo- 
hedral phase (i-FeF^; SG R3C) is antiferromagnetically ordered below T/v = 
363K while the amorphous varieties display a speromagnetic behavior below 
Tp 30-40K p ^lTll[T2] . The presence of fluorine in those materials establishes 
superexchange magnetic interactions between magnetic ions. The strength 
of the coupling is linearly dependent on cos^(^), where 9 denotes the superex- 
change angle Fe — F — Fe while the sign is given by the orbital conflguration 
of the FeFFe bond: for example, 180° {Fe^~^ — F^ — Fe^^) is characterized 
by a strong antiferromagnetic coupling according to Kanamori-Goodenough 
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rules [131 E]. 

The magnetic properties of -FeFs have been intensively studied during the 
past two decades and can be described in terms of iron ring statistics with 
induced magnetic frustration fl\ IT5]. 

• the rhombohedral phase i-FeF^, which is a pseudo-cubic packing of 
weakly tilted octahedral FcFq units, is not frustrated; 

• the hexagonal tungsten bronze WTB-FeF^ which results from the su- 
perposition of magnetically frustrated planes composed of hexagonal 
and triangular octahedral units cycles ; 

• and the pyrochlore pjT-FeF^ phase which consists of a packing of octa- 
hedral units building corner-sharing tetrahedron is thus more frustrated 
than HTB-FeFg. 

Recent experimental and theoretical studies on the properties of FeF^ nanos- 
tructures obtained by mechanical milling have shown that, at the nanoscale, 
this material is composed of two parts: the grain which behaves as the crys- 
talline T-FeF^ phase and a disordered grain boundary which remains com- 
posed of corner-sharing octahedral units [16l [IT]. In terms of ring statis- 
tics, it has been numerically established that all odd rings are localized in 
the grain boundary of nanostructures, confirming the speromagnetic struc- 
ture evidenced from in-field Mossbauer spectrometry. The Monte Carlo- 
Metropolis method has been successfully applied to study the surface and 
finite-size effects in nanostructures [TBI [191 EQ]. Nevertheless, no theoretical 
studies on the magnetic properties of ionic nanostructures which focus on the 
role of grain boundary by Monte Carlo simulation have yet been reported, 
to our knowledge. Unlike micromagnetic or molecular field calculations, the 
Monte Carlo simulation can take into account the atomic structure of the 
lattice and the nature of superexchange interactions for ionic materials such 
as FeFs. 

The Monte Carlo-Metropolis simulated annealing technique is an effective 
approach in the study of a system with many degrees of freedom. Dur- 
ing such a simulation, random numbers are used to sample conformations 
of the system with correct thermodynamical probabilities. A typical Monte 
Carlo simulation consists of two steps: thermalization and sampling. Dur- 
ing thermalization, the system is led adiabatically to its thermodynamical 
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equilibrium. After the system reaches this equihbrium, properties of interest 
can be estimated by averaging over enough samples [211 122] . Besides, in the 
simulated annealing method, it is possible to find heuristically the lowest 
energy state of the system. 

We used, in this work, the standard Metropolis algorithm to simulate the 
magnetic behavior of iron fluoride nanostructure with an emphasis on the 
role of grain boundaries. 

2. Magnetic model 

Our sample model is chosen as a double grain boundary constrained in 
between grains which consists of pure perfectly crystalline FeF^ with a size 
ranging from 8 to 12 nanometers as experimentally measured in iron fluo- 
ride prepared by high energy ball milling. Those samples are obtained by 
Voronoi tessellation and structurally relaxed with a scheme based on a mod- 
ified Metropolis algorithm [l7j. It has been shown both experimentally and 
numerically in such a system that the grain boundaries are disordered, but 
still remain composed of corner- sharing octahedral units with odd and even 
iron rings, thus leading to magnetic frustration. The macroscopic thermody- 
namic properties, such as the temperature dependence of the magnetization, 
the specific heat and the magnetic susceptibility for our system, are obtained 
from a Heisenberg-type Hamiltonian, which in general contains several terms 
corresponding to different energy contribution: exchange, Zeeman, dipolar, 
anisotropy, magnetostatic, magnetoelastic and thermal energy. In our study, 
we only consider the exchange energy for the first part of the work and then 
we add Zeeman contribution for the study of hysteresis loops. The Hamilto- 
nian of the system is thus 

^1 

i=l ^ jeV 

V is the nearest neighborhood of site i , Jij are the exchange coupling con- 
stant. Si and Sj are spins corresponding to the i an j sites, /ib is the Bohr 
magneton, gfj is the gyromagnetic ratio and H^xt is the external magnetic 
field. 

The first step of our work consisted in the determination of the coupling 
constant Jij which is related to the superexchange angle by the relation 

J(e) = Ji8o° cos^ e + Jgoo sin^ 9 (2) 
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which can be written in the form 

J{9) = JgQo + (Ji80° - ^90°) COS^6' (3) 

where JgQo and Ji8o° are respectively the couphng constant corresponding to 
the superexchange angle of 90° and 180°. 

From previous studies several pieces of informations are available for r-FeF^ 



The Neel temperature is approximatively 363K [15j; 

The blank angle (angle for which the constant coupling is zero) is 115° 



• The superexchange angle in this phase is 153.15° 

Taking into account the fact that each iron atom has six F neighbors and 
that octahedral units are regular in pure perfectly crystalline i-FeF^, we can 
use the relation giving the Neel temperature as a function of the coupling 
constant to determine the value of J^^ in units of the Boltzmann constant 

Tn = 1-45^^ (4) 
Qkb 

where Z is the coordination number and S the spin. With Z = 6 and S = 1 
Jq^ = -250.35 K 

By combining relation [3] for the rhombohedral superexchange angle (Or) and 
the blank angle (Ob), we obtain the system 



J{6r) = ^90° (1 - COS^ Or) + ^180° cos^ Or 
= J90°(l — COS^ 9b) + ^180° cos^ 9b 



(5) 



The resolution of this system gives 



_ Jeji cos^ 9b 

-'90° — — ^ (,DJ 

cos^ 9b — cos"^ 9r 

and 

_ J{9r){cos'9b-1) 

-'180° — 5-^5 ^Tn \') 

COS^ 9b — COS^ 9r 

which gives Jqqo = 72A2K and Jisqo = — 333.07/^. 

We can then deduce the relation giving the coupling constant as a function 
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of the superexchange angle. This relation is represented in figure [T} A com- 
parison with the distribution of the superexchange angle in nanostructured 
FeF^ obtained from previous studies [IT] implies that all interactions in the 
GB remain antiferromagnetic. This is observed in figure [T] by the large peak 
centered around 160°, all angles being greater than 140°. It is thus concluded 
that the frustration in the GB does only originate from the cationic topology. 

J{e) = 72.42 - 405.49 cos^ 9 (8) 
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Figure 1: (a): Coupling constant as a function of the superexchange angle for iron fluo- 
ride, (b): Angular distribution function for the superexchange angle. GB is characterized 
by a large peak around 160. 

In the first step of the simulations we only considered the exchange contri- 
bution to the magnetic Hamiltonian of the system. We have then neglected 
other phenomenological contributions such as dipolar, Zeeman or anisotropy 
terms, in order to clarify the infiuence of the exchange coupling on the mag- 
netic behavior. All samples, with sizes ranging from 8 to 12 nanometers and 
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different orientations of grains (tilted, twisted etc) were first structurally re- 
laxed. The energy given in equationjl] is minimized by means of the Monte 
Carlo/Metropolis simulated annealing procedure. During the simulation, 
Monte Carlo steps are applied on individual magnetic moments ("spins") 
with a random walk while solid angles are uniformly distributed over Air. 
Starting with a random spin configuration at a temperature Tg^gj.^ much 
higher than T/v (600K), the energy is minimized using the simulated anneal- 
ing scheme with a decreasing power law for temperature T = ^^startO-^'''* 
where i is the step number, while the thermodynamic quantities, such as 
magnetization, susceptibility and specific heat, can be derived as a function 
of the temperature. In each simulation, the final temperature is lower than 
IK and the number of Monte Carlo steps per spin (MCS) is 10^. 



3. Results 

We first determine the magnetization, susceptibility and specific heat of 
the system as a function of temperature. Figure |2] represents the specific 
heat and one can observe a transition around 360 K which is the Neel tem- 
perature of r-FeF3. The magnetization of the system is uniformly zero due 
to antiferromagnetic behavior. The magnetic configuration after annealing 
is composed of two antiferromagnetic grains separated by a transition zone 
with disordered spin structure as illustrated on figure [2j Each change in 
simulation conditions (number of MCS, orientation of initial grains, final 
temperature, etc) leads to a different orientation of magnetic moments of 
grains and then a different coupling of GB moment. 

We can conclude that the system is not in the magnetic equilibrium state. 
This can be easily understood by considering the fact that the GB is mag- 
netically frustrated. The present single spin flip we used is not appropriate 
for finding the magnetic equilibrium state in a reasonable time of simulation. 
New conditions have to be applied to reach this magnetic equilibrium state. 
A constraint can be added by considering the fact that there is one orien- 
tation of the relative global moment of each grain for which the magnetic 
energy is minimal. We used fixed boundary conditions to constrain magnetic 
moments of the two grains to form an angle 6 and then relaxed the interface. 
Figure |3] represents the relative energy of the system for an angle varying 
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Figure 2: (a): specific heat of the system as a function of temperature (Neel temperature 
is around 360K).(h): Magnetic configuration after annealing: each change in simulation 
conditions would give a different orientation of the total moment of the grain. 

from 0° to 360°. We can clearly distinguish two magnetic states with minimal 
energy around 90° and 270°. 

The same behavior of the magnetic energy was obtained for different grain 
orientations, but the angle corresponding to the minimal energy was not the 
same. Figure [3] compares the same results obtained with non-tilted, tilted 
and twisted samples : it clearly appears that the angle corresponding to the 
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Figure 3: (a):Magnetic energy as a function affixed angle boundary condition. There are 
two magnetic states of lower energy around 90° and 270°. (b): The same behavior is 
obtained for tilted, twisted and non-tilted system but with different angles. 

minimum energy depends on the initial relative crystallographic orientation 
of grains. It has recently been shown [17] that the only difference between 
those samples is the rate of frustration in the GB, suggesting that the angle 
corresponding to the minimum energy depends on the rate of frustration. 
To verify this dependence, we first considered a simple model consisting of a 
cubic cluster with all interactions being antiferromagnetic, then we changed 
some antiferromagnetic in ferromagnetic interactions at the interface to cre- 
ate locally a bond frustration. Frustration rate can be thus defined as 
T = NJf/{NJaf + NJp), where NJp and NJaf respectively represent the 
number of ferromagnetic and antiferromagnetic bonds. For each rate, the 
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system is relaxed and the angle between the moments of both sides of inter- 
face is measured: a certain rate of frustration gives rise to a non collinear 
orientation of the moments as evidenced on figure |4} 

Therefore, we can assume that the non collinear orientation observed in 
ferric fluoride nanostructures is due to the magnetic frustration of the grain 
boundary. We conclude that the angle between the total momenta of the two 
grains depends on the rate of frustration, which suggests the possibility to 
tune the final orientation of the grains. To achieve the study of this system. 




Figure 4: Coupling angle between grains as a function of frustration rate. Non collinear 
coupling appears for a rate comprised between 0. 2 and 0. 8. 

we have performed a simulation of hysteresis loops. It is well established 
that crystalline i-FeF^ is antiferromagnetic and that the GB behaves as a 
speromagnet giving rise to zero magnetization : therefore we do not expect 
any remanence on the hysteresis loops. The goal of the study is to analyze 
the influence of disorder in the GB on the spin rotation under a magnetic 
field. 

The hysteresis loops have been computed starting from a saturated state 
achieved after the application of a high enough field (2700 T) along the z- 
axis (orthogonal to the GB plane) and decreasing the field in constant steps 
6h = lOT, during which the magnetization was averaged over 5000 MC steps 
after thermalization. The temperature was sufficiently low (10~^ K) to avoid 
thermal fluctuations. The result for a sample with a 12 nm size is presented 
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in figure |5} As expected, there is no remanent magnetization in hysteresis 
loops, but the spin rotation from a configuration parallel to the magnetic field 
to an antiferromagnetic configuration is not uniformly done: the moments 
far away from the interface are decoupling first due to their strong coupling 
compared to those of the GB. 




Figure 5: Hysteresis loops and snapshots of the magnetic configuration during spin rotation 
for a 12nm size nana structure. As expected, there is no remanent magnetization hut a non 
uniform spin rotation in the area in inset. As the strength of the magnetic coupling far 
from the GB is greater than the one in the GB, magnetic moments in the GB are still in 
the direction of applied field after the spins far from the GB have rotated. 



4. Conclusion 

We presented a Monte Carlo study of nanostructured ferric fluoride by 
coupling a previously established structural atomistic modeling to a magnetic 
modeling. As experimentally evidenced, at the nanoscale, this compound is 
composed of two parts: A perfect antiferromagnetic configuration in the 
grain and a disordered magnetic configuration (speromagnetic like) at the 
grain boundary. The present numeric study has shown that the magnetic 
configuration at the grain boundaries is governed by the cationic topology 
which results from the disorientation of crystalline neighboring grains. Conse- 
quently, the lowest magnetic energy of the total system is strongly correlated 
to the rate of magnetic frustration in the GB, opening the possibility to tune 
the final orientation of the magnetic moments of the grains by modulating 
the frustration rate. 

From the careful simulation of hysteresis loops, the magnetization rever- 
sal is not uniform under a magnetic field ; in addition it is dependent on 
the strength of the magnetic coupling which is related to the superexchange 
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angle. It is finally important to emphasize that such an approach to numer- 
ically model structural and magnetic structures can easily be extended to 
ionic nanostructures, particularly to multilayers composed of different oxide 
species in order to follow the dependence of their properties on their topo- 
logical ordering, the thickness of the interfaces and the mismatch between 
layers. 
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